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Abstract — This paper considers broadcast ctiannels witii L 
antennas at tlie base station and m single-antenna users, wliere 
each user has perfect channel Itnowledge and the base station 
obtains channel information through a finite rate feedback. The 
key observation of this paper is that the optimal number of on- 
users (users turned on), say s, is a function of signal-to-noise ratio 
(SNR) and other system parameters. Towards this observation, 
we use asymptotic analysis to guide the design of feedback and 
transmission strategies. As L, m and the feedback rates approach 
infinity linearly, we derive the asymptotic optimal feedback 
strategy and a realistic criterion to decide which users should 
be turned on. Define the corresponding asymptotic throughput 
per antenna as the spatial efficiency. It is a function of the number 
of on-users s, and therefore, s should be appropriately chosen. 
Based on the above asymptotic results, we also develop a scheme 
for a system with finite many antennas and users. Compared with 
other works where s is presumed constant, our scheme achieves a 
significant gain by choosing the appropriate s. Furthermore, our 
analysis and scheme is valid for heterogeneous systems where 
different users may have different path loss coefficients and 
feedback rates. 
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I. Introduction 

It is well known that multiple antennas can improve the 
spectral efficiency. This paper considers broadcast channels 
with L antennas at the base station and m single-antenna users. 
To achieve the full benefit, perfect channel state information 
(CSI) is required at both receiver and transmitter. Perfect CSI 
at the receiver can be obtained by estimation from the received 
signal. However, if CSI at the transmitter (CSIT) is obtained 
from feedback, perfect CSIT requires an infinite feedback rate. 
As this is not feasible in practice, it is important to analyze 
the effect of finite rate feedback and design efficient strategy 
accordingly. 

The feedback models for broadcast channels are described 
as follows. To save feedback rate on power control, we assume 
a power on/off strategy where each user is either turned on 
with a constant power or turned off, and that the number of on- 
users (the users turned on) is a constant, say s, independent of 
the channel realization. For any given channel realization, the 
users quantize their channel states into finite bits and feedback 
the corresponding indices to the base station. After receiving 
the feedback from users, the base station decides which users 
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should be turned on and then forms beamforming vectors for 
transmission. 

Broadcast channels with feedback have been studied in [1], 
[2]. Ideally, if the base station has the perfect CSI, zero-forcing 
transmission avoids interference among users. However, with 
only finite rate feedback on CSI, the base station does not 
know the perfect channel state information and therefore in- 
terference from other users is inevitable. The interference gets 
so strong at high signal-to-noise ratio (SNR) regions that the 
system throughput is upper bounded by a constant even when 
SNR approach infinity. This phenomenon is called interference 
domination and was reported on in [1], [2]. The analysis is 
based on the assumption that the number of on-users s always 
equals to the number of antennas at the base station L (L < m 
is typically assumed). To limit the interference to a desired 
level, Sharif and Hassibi let m grow exponentially with L such 
that there are L near orthogonal users with high probability [1]. 
In both scenarios, a homogeneous system is assumed where 
all the users share the same path loss coefficient and feedback 
resource. 

Different from the above approaches, this paper studies het- 
erogeneous broadcast systems, where different users may have 
different path loss coefficients and feedback rates. Further- 
more, different from [1], we focus on systems with a relatively 
small number of users. Note that a cooperative communication 
network can often be viewed as a composition of multi-access 
and broadcast sub-systems with a small number of users. 
Research on broadcast systems of small size provides insights 
into cooperative communications. 

For such systems, we solve the interference domination 
problem by choosing the appropriate number of on-users s. 
The reason that random beams construction in [1] fails in our 
small size systems is elaborated in Theorem |2] 

Our solution is based on the asymptotic analysis where 
L, TO, s and the feedback rates approach infinity linearly with 
constant ratios among them. This type of asymptotics is 
applied to systems of small size. The main asymptotic results 
are: 

> It is asymptotically optimal to quantize the channel direc- 
tions only and ignore the channel magnitude information. 
The asymptotically optimal feedback function and code- 
book are derived accordingly. 

« A realistic on/off criterion is proposed to decide which 
users should be turned on. 

• The corresponding throughput per antenna converges to a 



2 



constant, defined as the spatial efficiency. It is a function 
of the normalized number of on-users s — j-. Further, 
there exists a unique s G (0, 1) to maximize the the spatial 
efficiency. 

We develop a scheme to choose the appropriate s for 
systems with finite L and m. Simulations show that the 
gain achieved by choosing s is significant compared with 
the strategies where s = L. In addition, our scheme has the 
following advantages. 

• It is valid for heterogeneous systems. 

• The set of on-users is independent of the channel real- 
ization. As a result, computation complexity is low since 
we do not have to perform a user selection computation 
every fading block. 

• Only on-users need to feedback CSI, which saves a large 
amount of feedback resource. 

This paper is organized as follows. The system model 
is introduced in Section HI] Then Section |III] performs the 
asymptotic analysis obtaining insights into system design, and 
quantifies the spatial efficiency. Based on the asymptotic re- 
sults, a practical scheme is developed in Section|IV]for systems 
with finite many antennas and users. Finally, conclusions are 
summarized in Section FVl 

II. System Model 

Consider a broadcast channel with L antennas at the base 
station and m single-antenna users. Assume that the base 
station employs zero forcing transmitter Let 7^ G M\K~ 
(1 < i < m) be the path loss coefficient for user i. Then 
the signal model for user i is 



^ ^ ^ is the zero-forcing 
C is the source signal 



where G C is the received signal for user i, h; G 
the channel state vector for user, G i 
beamforming vector for user j, Xj G 
for the user j and Wi G C is the complex Gaussian noise 
with zero mean and unit variance CAf {0, 1). Here, we assume 
that q^Qj — 1 and the Rayleigh block fading channel model: 
the entries of hj are independent and identically distributed 
(i.i.d.) CAf {0,1). Without loss of generality, we assume that 
L < m; if L > m, adding L — m users with 7^ = yields an 
equivalent system with L' — m. 

For the above broadcast system, it is natural to assume 



total power constraint J2iLi^ 



< p. Further, for 



implementation simplicity, we assume a power on/off strategy 
with a constant number of on-users as follows. 

Al) Power on/off strategy: a source Xi is either turned 
on with a constant power Pon or turned off. It is 
motivated by the fact that this strategy is near optimal 
for single user MIMO system [3]. 
A2) A constant number of on-users: we assume that the 
number of on-users s(l < s < m)isa constant 
independent of the specific channel realizations. With 
this assumption, Pon = f • Here, s is a function of 
p, 7i and feedback rate. This assumption is different 



from the one in [1], [2], where s — L always (L < m 

is assumed there). 
The finite rate feedback model is then described as follows. 
Assume that both base station and user i knows 7J1] but 
only user i knows the channel state realization perfectly. 
For given channel realizations hi • • • h^, user i quantizes his 
channel hj into Ri bits and then feeds the corresponding index 
to the base station. Formally, let Bi — |h ^ C^^^j with 

\Bi\ ~ 2^' be a channel state codebook for user i. Then the 
quantization function is given by 



q : C 



Lxl 



In Section ITlI-AI and Illl-Bl we will show how to design q and 
B respectively. 

After receiving feedback information from users, the base 
station decides which s users should be turned on and forms 
zero-forcing beamforming vectors for them. Let Aon be the set 
of s on-users. The zero-forcing beamforming vectors q/s i G 
Aon is calculated as follows. Let V^- be the plane generated by 
|hj : j G j4on\ {*}|- Let Vi be the orthogonal complement 
of Vt and t be the dimensions of P,. Define T, G C^""* the 
matrix whose columns are orthonormal and span the plane P. 
Then q^ is the unitary projection of h^on 

T T^h 

q« := r:— • (1) 



T,T|h, 

Here, if s = 1 and Aon — is a L x L unitary matrix 

and qi = hi/ h^ 

III. Asymptotic Analysis 

In order to obtain insights into system design, this section 
performs asymptotic analysis by letting L,m,R/s 00 
linearly. The quantization function q and asymptotically op- 
timal codebook Bi are derived in Section IIII-AI and IIII-BI 
respectively. Then Section IIII-CI develops a realistic on/off 
criterion to decide which users should be turned on. Finally 
Section HlI-DI computes the corresponding spatial efficiency. 

A. Design of Quantization Function 

Generally speaking, full information of h^ contains the 
direction information -.^ hj/||hi|| and the magnitude 
information ||hi||. In our Rayleigh fading channel model, it 
is well known that v; and ||hi|| are independent. Intuitively, 
joint quantization of Vj and ||hj|| is preferred. 

Interestingly, Theorem [T] implies that there is no need to 
quantize the channel magnitudes. Indeed, as L,m 00 
linearly, all users' channel magnitudes concentrate on a single 
value with probability one. 

Theorem 1: For Ve > 0, as P, m ^ 00 with ^ fh ^ K+, 



Pr 



max — II hj 

l<i<7n L 



> 1 



0, 



'There are many ways in which the base station obtains 7^. A simple 
example could be that the base station measures the feedback signal strength. 
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and 



\l<i<m Li J 

The proof of Theorem [1] is omitted due to the space 
Hmitation. An important fact behind the proof is that whether 
the users' channel magnitudes concentrate or not depends on 
the relationship between L and m: this concentration happens 
in our asymptotic region where L and m are of the same order 

To fully understand Theorem [1] it is important to realize 
that the Law of Large Numbers does not imply that all users' 
channel magnitudes will concentrate. According to the Law 
of Large Numbers, ||hi|| 1 almost surely for any given 
i. However, if m approaches infinity exponentially with L, 
there are certain number of users whose channel magnitudes 
are larger than others', and therefore it may be still beneficial 
to quantize and feedback channel magnitude information. This 
phenomenon is illustrated by the following example. 

Example 1: (A case where magnitude information is bene- 
ficial) Consider a broadcast channel with 71 = • • • = 7„i = 1. 



As L, m ^ cx) with log (m) / L m! 
e > 0, (5i > and 82 > such that 



€ ffi."*", there exists an 



L 



log 



1 
L 



and 



^log 



^ llh. 



5i. 



with probability one. Note that there are a set of users whose 
channel magnitudes are e-larger than another set of users. It 
may be worth to let the base station know which users have 
stronger channels. 

Theorem [U implies that it is sufficient to quantize the chan- 
nel direction information only and omit the channel magnitude 
information. For this quantization, the codebook is given by 
Si = {p e C^^i : IIpII = 1} with \B^\ = 2^-. We use the 
following quantization function 

q : C^^i 



hi I— > Pi = arg max 
pes, 

where Vi is the channel direction vector 



v|p 



(2) 



B. Asymptotically Optimal Codebooks 

Given the quantization function (|2]i, the distortion of a 
given codebook Bi is the average chordal distance between 
the actual and quantized channel directions corresponding to 
the codebook Bi and defined as 



D {Bi) := 1 - Eh. 



The following lemma bounds the minimum achievable dis- 
tortion for a given codebook rate (usually called the distortion 
rate function). 

Lemma 1: Define D* (R) = inf D (B). Then 

B: |e|<2« 







2" 








max 






pes. 







1 



2'— (1 + 0(1)) < D* {R) 



< 



L - 1 



2-— (1 + 0(1)) 



(3) 



and as L and R approach infinity with 



r G 



lim D* (R) = 2-^. 

The following Lemma shows that a random codebook is 
asymptotically optimal with probability one. 

Lemma 2: Let Siand be a random codebook where the 
vectors p e brand's are independently generated from the 
isotropic distribution. Let R = log|Srand|- As L, R ^ 00 
with ^ ^ f e IR+, for Ve > 0, 



> 2-^ + e} = 0. 



lim Pr {Brand : D {B^- 

a 

The proofs of Lemma [T] and |2] are given in our paper [4]. 

Due to the asymptotic optimality of random codebooks, we 
assume that the codebooks Si's z = 1, • • • ,m are independent 
and randomly constructed throughout this paper. 



C. On/ojf Criterion 

After receiving feedback from users, the base station should 
decide which s users should be turned on. 

Ideally, for given channel realizations hi, - - , h„i, the 
optimal set of on users A^^^ should be chosen to maximize 
the instantaneous mutual information. Note that the base 
station only knows the quantized version of channel states 
Pir ' iPm- It can only estimate the instantaneous mutual 
information through Pi's. The set is given by 



^on = arg max ^ log (1- 

Ao„: \A^„\=s ^^^^^ 



7if plqi 



(4) 



However, finding ^4*,^ requires exhaustive search, whose com- 
plexity exponentially increases with to. 

The random orthonormal beams construction method in [ 1 ] 
does not work for our asymptotically large system either In 
[1], the base station randomly constructs L orthonormal beams 
bi, - - - ,bi, finds the users with highest signal-to-noise-plus- 
interference ratios (SINRs) through feedback from users, and 
then transmits to these selected users. There, the SINR cal- 
culation for user i is related to the quantity max hjbi, . 

l<k<L 

However, Theorem |2] below shows that in the asymptotic 
region where L and to are of the same order, all users' 
channels are near orthogonal to all of the L orthonormal beams 
bi's. Therefore, all users' maximum SINRs (maximum over 
L given orthonormal beams) approach zero uniformly with 
probability one. The method in [1] fails in our asymptotically 
large system. 

Theorem 2: Given Ve > and any L orthonormal beams 
bfe G C^^i 1 < fc < L, as L, TO ^ 00 linearly with ^ 
fh e K+, 



lim Pr 



1 

max — 

l<i<m, l<k<L L 



> e 



0. 



The proof is omitted due to the space limitation. 
In this paper, we take another approach where the on/off 
decision is independent of channel directions. We start with 
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the throughput analysis for a specific on-user i E Aon- Note 
that 

The signal power and interference power for user i are given 



by 



and 



Slg,2 



(5) 



(6) 



jeAo„\{i} 

respectively. Note that the influence of the users in Aon\ {*} on 
user i only occurs through their directions q^'s j £ Ac,n\ {«}■ 
If the choice of on-users Aon is independent of their channel 
directions, then and q^'s are independent. In this case, -Psig,i 
and Pint,i can be quantified as L, m, s, R/s — > oo. The result 
is given in the following proposition. 

Proposition 1: Let \Aon\ — s and L,m, s, R/s oo with 
^ ^ TO, I; — > s and — > fi. Assume that v^'s i e Aon are 
independent. Then for Vi G ^on. 



P • 

Slg,2 



^7.(1 

s ^ 



O(i-s-) 



p 



int.i 



and therefore 



X,; := log ( 1 
^ log ( 1 



P 



1 + Pint,i 
1 - S 



with probability one, where 







(7) 



1 + pji2-^^ 

Remark 1: This proposition may not be true if Vj's (j G 
^on\ {*}) are not independent of v^. Indeed, for example, if 
other users are chosen such that their channel directions are 
as orthogonal to user i as possible, the interference to user i is 
less than that achieved by our choice where channel directions 
are not taken into consideration. This claim is verified by the 
fact that 3e > such that 



mm 

Ao„\{i:} 



je^o„\{n 



< E 

jeAo„.„„d\{i} 



hlq, 



with probability one as L,m,s oo linearly, where 

^on.rand\ {*} dcnotcs a random choice of Aon\ {*}■ 

Remark 2: Proposition [T] shows that the user i's asymptotic 
throughput is a constant independent of the specific channel 
realization hi with probability one. 

Based on Proposition [T] we select the set of s on-users Aon 
such that I Aon I — s and 



Ao 



{i : rii> r/j for \fj ^ Aon} ■ 



(8) 



if there are multiple candidates, we randomly choose one 
of them. It is the asymptotically optimal on/off selection if 
the on/off decision is independent of the channel direction 



information. The difference between the throughput achieved 
by optimal on/off criterion in (|8]l and the proposed one in (|4]i 
remains unknown. 



D. The Spatial Efficiency 

We define the spatial efficiency (bits/sec/Hz/antenna) as 

X(s):= lim 

where L, to, s, R/s ^ oo in the same way as before, X^^^ is 
the average throughput per antenna given by 



and Aon, Psig,i and Pint,i are defined in (|8]l, ^ and Q 
respectively. 

We shall quantify T (s) for a given s. Define the empirical 
distribution of rji as 

M(")(r/<x) ■.= -\{nr. m<^}\, 

' m 

and assume that '■— lim/x^™'' exists weakly as 

L, TO, Pi's oo. In order to cope with /i,,'s with mass points, 
define 



for Vcc S R, where / is a integrable function with respect to 
jirj. Then T is) is computed in the following theorem. 



Theorem 3: Let L,m, s, Ri's oo with ^ 
fi. Define 



TO, 



and ^ 



ris := sup <^r/ : m I d^iri > s 



Then as s ^ (0, 1), I (s) = 0. If s e (0, 1), 
X{s) = fh I log ( 1 + ?; — — ) d/^^ 




dur^ log 1 + 77. 



1 



(9) 



We are also interested in finding the optimal s to maximize 
T{s). Unfortunately, I (s) is not a concave function of s 
in general. Furthermore, the measure /i^ may contain mass 
points. The optimization of I (s) is therefore a non-convex 
and non-smooth optimization problem. The following theorem 
provides a criterion to find the optimal X(s). 

Theorem 4: T. (s) is maximized at a unique s* e (0, 1) such 
that 



e 



lim mi — , 

As^O As 

X{s*)-X{s* ~ As) 



lim sup 



As 



(10) 



The proof is omitted due to the space limitation. The X (s*) 
is the maximum achievable spatial efficiency for the proposed 
power on/off strategy. 
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IV. Finite Dimensional System Design 

Based on the asymptotic results in Theorem [3]|4l we now 
propose a scheme for systems with finite L and m. 

A. Throughput Estimation for Finite Dimensional Systems 

While asymptotic analysis provide many insights, we do 
not apply asymptotic results directly for a finite dimensional 
system. The reason is that in asymptotic analysis 7; ^ 
while J- cannot be ignored for a system with small L. In 
the following, we first calculate the main order term of the 
throughput for user i e Aon and then explain the difference 
between asymptotic analysis and finite dimensional system 
analysis explicitly. 

To obtain the main order term, proceed as follows. Note 
that the throughput for user i e Aon (|Aon| — s) is 



E 



log 1 



= log 1 



1 + Pint,, 
1 + E [Fint.d 



E 



E 



log 

log 



1 ~^ PsiK.i ~^ P'li 



1 + E [P,i 

1 + Pint. 



E [Pi, 



where P,; 



,l+E[Pint,^] 

and Pint,i are defined in (|5]l and The fol- 
lowing theorem calculates E [Psig,i] and E [Pint,i] for finite 
dimensional systems. 

Theorem 5: Let S,;'s be randomly constructed and Di = 
Eg; [D (Bi)] for all 1 < i < m. For randomly chosen Aon 
and i £ Aon, 

- 1' 



E[Psi 



and 



if 1 < s < L 

L ' 



1 



L 



1 



P (i - 1) 

(fl) 



E[Pn 



L 



1 



(12) 



if s > P, E [P,ig,,] = 0. 

The calculation of E [Psig.i] and E [Pint.i] relies on quantifi- 
cation of Di. In general, it is difficult to compute P,; precisely. 
Note that the upper bound in (O is derived by evaluating 
the average performance of random codebooks (see [4] for 
details). We use its main order term to estimate Di. 



r 



P, 



1 

L-l 



Define 



log 1 



E[P.i 



1 + E [Pint A 

It can be verified from Proposition [T] that P; = 



(13) 

0(1) 



Then the difference between asymptotic analysis and finite 
dimensional systems analysis is clear. In the limit, s 
and 1^ f,. However, for finite dimensional systems, 
simply substituting these asymptotic values into ( II 1H13| | di- 
rectly introduces unpleasant error, especially when L is small. 



Therefore, to estimate Ij (Vi G Aon) for finite dimensional 
systems, we have to rely on ( I11H13I I. 

B. A Scheme for Finite Dimensional Systems 

Given system parameters P, to, ji's and P^'s, a practical 
scheme needs to calculate the appropriate s and Aon. This 
process is described in the following. 

For a given s, the set of Aon is decided as follows: we 
first calculate Pmain,i, ■ • ■ j^main.m according to ( HjI ) and then 
choose the s users with the largest Pmain.i's to turn on; if 
there exists any ambiguity, random selection is employed to 

resolve it. For example, if Pmain,l > Pmain,2 > • ■ ■ > Pmain.m, 

the user 1,2,- •• ,s are turned on. If Pmain.i = 2niain,2 = 
• • • = Pmain.m, the s on-users are randomly selected from all 
the m users. Note again, Aon is independent of the channel 
realization. 

The appropriate s is chosen as follows. Let 



max y Pmain,i- 



Here, note that Pmain.i is a function of s. For a given broadcast 
system, we choose the number of on-users to be 

■^main ^ ^"^^ max Pniain (s) • 

l<s<L 

Although the above procedure involves exhaustive search, 
the corresponding complexity is actually low. First, the calcu- 
lations are independent of instantaneous channel realizations. 
Only system parameters P, m, 7i's and P^'s are needed. 
Provided that 7i's change slowly, the base station does not 
need to recalculate s*nain and Aon frequently. Second, Ri — 
Rj in most systems. For such systems and a given s, the s 
on-users are just simply the users with the largest ji's. 

After calculating s*,ain and Aon, the base station broadcast 
Aon to all the users. For each fading block, the system works 
as follows. 

• At the beginning of each fading block, the base station 
broadcasts a single channel training sequence to help all 
the users estimate their channel states h/s. 

• After estimating their h^'s, the on-users quantize h^'s into 
Pi's according to (|2]l and feed the corresponding indices 
to the base station. 

• The base station then calculates the transmit beamforming 
vectors q^'s according to ([T]i, and then transmits q^X^'s. 

Remark 3 (Fairness Scheduling): For systems with 7^ 7^ 
7j or Ri Rj, there may be some users always turned 
off according to the above scheme. Fairness scheduling is 
therefore needed to ensure fairness of the system. There are 
many ways to perform fairness scheduling. Since fairness is 
not the primary concern of this paper, we only give an example 
as follows. Given ni users, the base station calculates the 
corresponding s*j^jj^ and Aon, and then turns on the users in 
Aon for the first fading block. At the second fading block, the 
base station considers the users who have not been turned on 
{I,-- - ,m}\Aon. It calculates the corresponding s'^^^^^ and 
Aon, and then turns on the users in the new Aon. Proceed this 
process until all users have been turned on once. Then start a 
new scheduling cycle. 
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L=4, m=4, Rfb=6 Bits/Channel Realization 

2.5 r 

Simulation : l(s„„. ) 
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(a) Rfb = 6 Bits/Channel Realization 



L=4, m=4, Rfb=12 Bits/Ctiannel Realization 
Simulation : l(s , ) 
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SNR (dB) 



(b) ijfb = 12 Bits/Channel Realization 
Fig. 1. Total Throughput for Zero Forcing Beamforming 



p and i?fb). The throughput of schemes with presumed s 
is presented in dotted Hnes. From the simulation resuhs, the 
throughput achieved by choosing appropriate s is always better 
than or equals to that with presumed s. Specifically, compared 
to the scheme in [2] where s = L — A always, our scheme 
achieves a significant gain at high SNR by turning off some 
users. 

V. Conclusion 

This paper considers heterogeneous broadcast systems with 
a relatively small number of users. Asymptotic analysis where 
L,m,s,Ri oo linearly is employed to get insight into 
system design. Based on the asymptotic analysis, we derive the 
asymptotically optimal feedback strategy, propose a realistic 
on/off criterion, and quantify the spacial efficiency. The key 
observation is that the number of on-users should be appro- 
priately chosen as a function of system parameters. Finally, a 
practical scheme is developed for finite dimensional systems. 
Simulations show that this scheme achieves a significant gain 
compared with previously studied schemes with presumed 
number of on-users. 
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C. Simulation Results 

Fig. [U gives the simulation results for the proposed scheme 
using zero-forcing. In the simulations, L = rn — 4. For 
simplicity, we assume that 71 = 72 = • • • = 7m = 1 and 
Ri = R2 = ■ ■ ■ — Rm = Rih- With these assumptions, 
the s on-users can be randomly chosen from all the m users. 
Without loss of generality, we assume that Aon = {1, • • • ,s}. 
Let I (s) = J2i£A -^i- ffl '■^^ solid lines are the simu- 

lations of T(s*-„) while the dashed lines are the theoretical 

N Illdlll ' 

calculation of Xmain (Si'^iain)- The simulation results show that 
the optimal s is a function of p and i?fb. For example, s — 1 
is optimal when p e [15, 20]dB and i?fb = 6 bits, while s = 3 
is optimal for the same SNR region as i?fb increases to 12 
bits. The reason behind it is that the interference introduced 
by finite rate quantization is larger when i?fb is smaller: when 
i?fb is small, the base station needs to turn off some users to 
avoid strong interference as SNR gets very large. 

We also compare our scheme with the schemes where the 
number of on-users is a presumed constant (independent of 



